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Abstract 



We present here a modified hydrodynamical model which can be used to 
investigate in a easy way the collective electronic excitations of spherical 
atomic clusters. The model is then applied to the case of the giant-dipole 
plasmon mode of the free fullerene molecule. 



PACS numbers: 71.45. Gm, 36.40.+d 



Collective electronic excitations of small metallic particles or atomic clusters like 
fullerites may be and in fact are currently [1-6] of great interest with the goal to bridge 
the gap between our description of atom and solid. A part the polarization propagator 
response method |7j within RPA || || a good tool to investigate theoretically the 
plasmon-like collective behaviour of the nearly free or valence electrons in such systems 
is certainly the hydrodynamical model which date back to pioneering work by Bloch 



1C|. It was early used [ I I] eventually in refined version, as for instance in the 



investigation of photoabsorption of a Thomas- Fermi atom [ITSl . More recently such an 
approach was successfully used investigating the excitation of plasmons by energetic 
electrons in small spheres |13] and photon bombardment of a plasma slab JlIJ . Since 
then the model has been refined, fl5| and now leads to very sophisticated equations 
for the density fluctuations around the equilibrium distribution of a inhomogeneous 
electron gas |16| . In spite of that, this approach seems not so widely used, probably 
because the necessary equations are very cumbersome to handle unless the equilibrium 
density distribution can be considered constant. Thus I propose in this note a way to 
simplify these difficulties in the case of a problem with spherical symmetry, which will 
be presented within the more simple Bloch's original scheme. It consists on the idea 
to consider the equilibrium electronic density po( r ) which is isotropic for a problem 
with spherical symmetry, to assume constant values p^ within each shell of a given 
thickness A, by which the space around the origin can be subdivided. I will also 
consider the quantity above to be zero outside a certain distance R form the origin. 
The case of po(r) being continuous is of course recovered if the shell's thickness goes to 
zero and R, if necessary, to infinity. But all that amounts to present an approximate 
method to solve the continuous case, with no less difficulties than before if you like 
to approach more and more the limit. Instead, also in view of the validity of the 
approximation leading to po(r) both the set of values pg as well R and the Aj will 
be considered as parameters referring to a model equilibrium density. The net result 
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will be a kind of multi-shell or onion model. 

The starting point in the Bloch's approach is a Lagrangian £2- It comes as the 
second order (quadratic) contribution in the development around equilibrium of the 
full Lagrangian C for an irrotational fluid of electrons in which only the Coulomb in- 
teraction with the positive background, the Coulomb interaction among the electrons 
and the kinetic energy term are considered. The first variation of £ in the static case 
leads to the Thomas-Fermi equation. The equations of motion ( written in atomic 
units to be used in the following ) result to be: 

0(r, t) = n(r, t) + [ ^4 dr ' ™( r ' *) = V " [Po(*W(r, t)) (1) 

Po(r) J \r — r'\ 

where n(r,t) and 4>(r,t) are perturbations in the electronic density and velocity 
potential , respectively, in the electron gas. As already indicated po( r ) is the elec- 
tronic density in the undisturbed state so that po(r)'V(j)(r,t) represents at this or- 
der of development the current density of the fluctuations. C[po( r )] is the quantity 
7T 2 /3(3/7r) 2/3 p (y) 2/3 . Assuming: 

*('.«> = /(-){ = < 2 > 

we have : 

f , x C[p (r)] f g{r') , 

uf(r) = rrW+ /i 7\ dr 

Po(r) J \r — r'\ 

cog(r) = -V-[p (r)V/(r)]. (3) 

In the case that po( r ) assumes a constant value p and using both the Eq. (H) one see 
that g(r) satisfies the equation : 

A S (r) + (^2o) sW =0 (4) 

which is easy solvable after separation of the angular dependence. In the present case 
it is convenient to do so using the real, normalized spherical harmonics y^^Q) |~7| 
where s means e (even) or o (odd), in order to easy obtain f(r) and g(r) as real 



functions as it should be. Let F = (uj 2 — 4irp )/C[p ] and p = J\F\. The general 
solution of Eq. (|) is of the type: 

gl m (r) = [aAtOir) + bB l (pr)} 3^(0) (5) 

where a, b are numerical coefficient and Ai(pr), Bi(pr) indicate the couple ji{pr) and 
T]i{fir) of independent Bessel functions of order I when F is positive or their analytical 
extension in terms of hyperbolic sinus and cosinus otherwise. Consider now a isotropic 
equilibrium density distribution which is zero outside the sphere of radius R and let 
R be subdivided for simplicity in N equal intervals of length A = R/N. The sphere 
is thus subdivided in N shells, the i th among them having rj_i and as internal and 
external radius respectively with = iA (i = 1, N) and let the quantity p^ be the 
constant value of the equilibrium density within the i th shell. The solutions for g(r) 
will be: 

N 

9t m (r) = EK ( S MW) + B M\ y S m,iW X ( Hr) (6) 

i=l 

where x i r ) is the step function which is equal to one inside the i th shell and zero 
otherwise, and the meaning of the quantities indexed with i is obvious. The corre- 
sponding solution for the velocity potential is obtained from Eq. (H) using the first 
equation in Eq. (|]). It is required that both f(r) and g(r) be regular at the origin 
and as usual that the normal component of the fluctuation current density at the sur- 
face r = R be zero. Besides these boundary conditions also the matching conditions 
at the separating surfaces of the shells are required. For that we have at disposal 
the values of the fields and n and of their derivatives at such places. As po(r) is 
considered discontinuous I see no need to impose continuity conditions on the fluc- 
tuation density. So we can concentrate on requiring the continuity of the velocity 
potential and of the fluctuation density current po(r)'V<j)(r,t). The regularity of the 
fields at r = requires that the coefficient of the non regular Bessel function referring 
to the first shell around the origin be zero. The boundary condition at r = R and 
the 2N — 2 conditions at the separation walls then give rise to 2N — 1 homogeneous 
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linear equations for the remaining 2N — 1 unknown coefficients and we are left with 
an eigenvalue problem. Let u n be an eigenvalue and f n (r) g n {?) the couple of cor- 
responding solutions for f(r) and g(r). It is easy to show |18|] that thanks to the 



conditions imposed the integral of the product f m (r)g n (r) over the volume \r\ < R 
is zero unless m = n, that is the result which holds when po( r ) is continuous || [12 
is recovered. The solution will then be normalized with the following choice : 

fn(r) g m {r)dr = S nm u n . (7) 

I notice that a choice of an orthogonal set of harmonic functions for the angular part 
was already made, so the degeneracy connected with the spherical symmetry does 
not cause trouble in this respect. 

To illustrate the method outlined above I applied the onion model to the collective 
electronic excitations of fullerene, since this recently produced macromolecule is the 
subject of many investigation both on experimental and theoretical side. The Ceo is 
considered with good approximation |19| as a neutral system composed of nearly-free 
electrons with the background positive charge uniformly distributed over the surface of 
a sphere of radius rv As a starting point the Thomas-Fermi equilibrium distribution 
Potf{t) of the valence electrons, whose number was fixed to be 240 following the 
Ref. |p5fl , was computed. Then the following quantities were considered for each 



Poi = 77 I POTF{r)dVi 



shell: 

Vi 

where Vi is the volume of the i th shell. The required constant quantities pg were 
then obtained as p 0i time a normalization factor with the aim to ensure that the total 
number of electrons contained in the volume |r| < R be exactly 240. The eigenvalue 
problem was solved numerically for three different values of the fullerene radius r& 
in the only case for which 1 = 1. The interval of explored frequencies ranges from 7 
to 32 eV, in order to explore the discovered high-energy giant-dipole plasmon mode 
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of the free molecule [20|. We know that it exists a lower limit of distance given by 
1/Z where Z is the charge, for the applicability of Thomas- Fermi equation near a 
charged nucleus (in the present case the carbon nucleus screened by the two inner 
s-electrons) [pi]] , and therefore the detailed behaviour of the resulting equilibrium 
electronic density at shorter distance can be ignored. We also know that the ground 



i.e. 



state electronic distribution of a Thomas-Fermi atom has a tail going wrong |22 
it goes too much slowly to zero at infinity, and as a crude remedy we can cut off the 
obtained Potf{ t ) at some finite radius. Taking into account the considerations above 
the values A = 0.2 and R = 14.2 were adopted as sensible. The last figure is about 
two time the value of 6*60 radius i.e. 3.5 A ||23|l and so the equilibrium distribution 
ends at fairly the same radial distance from the positive charged surface. After the 
eigenvalues and the related eigenf unctions are known the next step is to build up the 
Hamiltonian which using these normal mode will look of course like a collection of non 
interacting harmonic oscillator hamiltonians. Let us restrict to the space of functions 
whose angular dependance is that of the usual Z = l,m = harmonic. Indicating 
with Sfc(r) and Wk(r) the k th solution with the above angular part for f(r) and g(r) 
respectively and expanding with such a basis the fields <ft and n it is easy to build up 
the contribution to the total Hamiltonian as coming from the subspace considered. 
It is also easy to write down the dipole interaction term when an external electric 
field E directed along the polar axis is considered. At last one see that the oscillator 
strength Xk of the k th dipole mode is Xk = y\ where 

Uk = j t cos 'dwkir) dr . (9) 

I recall that the functions Sk{r) and Wk(r) are normalized as indicated previously. 

The results of computations are collected in Table 1. One may see that the con- 
tribution of the dipole strength within the range 7-32 eV to the total figure of 240 
electrons required by the dipole sum rule is quite large. After the inspection of other 
solutions below this range and above up to 42 eV it appear that the sum rule is satis- 



fied in defect within about then percent. So it is very probably satisfied exactly taking 
into account larger frequencies and the computing errors as well. The contribution of 
the modes listed to the static polarizability of a free Cqq give a somewhat large figure, 



in agreement with the quantum-chemical result of reference [24|. The imaginary part 



of the dynamic polarizability for a free fullerene molecule was recently determined 



indirectly from photoionisation data on gaseous Ceo p0[ . The authors concluded for 
the presence of absorption peaks at about 7.8, 10, 13, 17 and 21 eV. The plasmon 
collective excitation modes are presumably strongly coupled to the single electron 
states and so the absorption line width can be very large. Referring to the middle 
column of the Table 1 we could thus suppose that the experimental maxima near 13 
and 17 eV are due to the two couples of resonances here found at 12.14, 14.37 eV 
and at 16.41, 18.17 eV respectively. We may then conclude that the prediction of our 
calculations in that of the three cases for which the fullerene radius is the nearest to 



the indicated value of 3.5 A|23] is in a good accordance with the experiments. I stress 
however that the aim was especially to present here the model and not to present 
very refined predictions. The conclusion above still seems to indicate that the method 
outlined does work well. It can be useful to look at other multipoles than the simple 
dipole and then investigate such things like the inelastic electron scattering and the 
behaviour of Van der Waals interaction at large distance to which large multipoles 
may contribute appreciably [^, . As a concluding remark I observe that of course 
the present method can be improved further. First such terms like exchange and 
correlation energies must be included in the Lagrangian of departure. And on the 
mathematical side it would be useful to correlate the shells thickness to the rapidity 
of variation of pg . 
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Table 1: The frequencies uj and the oscillator strengths x for three values of r b as 
result from the present calculation in the frequency range 7-32 eV. The cases where 
x < 0.25 are omitted. 



n = 


6.4 


n = 


6.6 


n = 


6.8 


uj (eV) 


X 


uj (eV) 


X 


uj (eV) 


X 


7.23 


9.76 


7.60 


11.26 


8.00 


13.11 


9.44 


10.44 


9.92 


12.46 


10.43 


15.25 


11.57 


14.34 


12.14 


18.00 


12.71 


21.98 


13.70 


18.03 


13.85 


1.41 


12.95 


1.39 


14.89 


0.69 


14.37 


22.32 


14.99 


33.62 


15.77 


25.19 


16.41 


36.73 


16.88 


37.40 


17.71 


35.55 


18.17 


29.75 


17.61 


5.79 


19.38 


28.30 


18.92 


7.76 


18.95 


32.81 


20.28 


4.66 


20.20 


27.08 


20.85 


7.51 


21.40 


23.54 


22.09 


6.36 


22.24 


5.78 


23.20 


6.84 


23.76 


2.82 


23.21 


4.59 


25.82 


7.83 


24.42 


6.32 


26.87 


4.94 


27.39 


1.00 


26.46 


0.36 


27.88 


0.37 


30.94 


4.38 


28.50 


1.50 


30.15 


0.25 


31.74 


0.58 


29.06 


3.60 


31.51 


4.18 
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